Quantum phases of disordered flatband lattice fractional quantum Hall systems 
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By numerical exact diagonalization techniques, we obtain the quantum phase diagram of the lat- 
tice fractional quantum Hall (FQH) systems in the presence of quenched disorder. By implementing 
an array of local potential traps representing the disorder, we show that the system undergoes a 
series of quantum phase transitions as the disorder and/or the interaction is tuned. As the strength 
of potential traps is increased, the FQH state turns into a compressible liquid, and then into a topo- 
logically trivial insulator. We use numerically calculated energy gap, quantum degeneracy, Chern 
number, entanglement spectrum, and fidelity to identify various quantum phases. The connection 
to continuum FQH effects is also discussed. 



PACS numbers: 71.10.Fd, 03.75. Ss, 05.30. Fk, 11.30.Er 

Introduction — As one of the most fascinating discover- 
ies in physics, the fractional quantum Hall efFect(FQHE) 
has been the central focus of intense research [l|-|3[ over 
the last three decades. Recently, the interest in this topo- 
logical state has been further enhanced by the theoretical 
proposal and the numerical discovery of a new kind of 
QHE in the absence of an external magnetic field, using 
lattices vifith nearly- flat topological bands [Jj^- Exten- 
sive work has established that these novel lattice FQH 
states are topologically identical to the ordinary contin- 
uum FQHE observed in 2D electron systems [7Hl5| . This 
recent theoretical realization of the lattice version of the 
well-established continuum QHE (and FQHE) has cre- 
ated great excitement in the community because it sub- 
stantially extends the range of possible topological phases 
which may exist in nature. 

Nearly three decades after its discovery, the topological 
nature of the FQHE is now reasonably well-understood. 
However, some key questions with considerable exper- 
imental implications remain open. Among these chal- 
lenges, perhaps the most important is understanding the 
effect of disorder. As a gapped state, a topologically non- 
trivial FQH phase is stable against weak disorder. As 
the strength of disorder increases, the FQH phase is sup- 
pressed and the system eventually becomes a trivial in- 
sulator. Although some transport signature has been de- 
tected in numerical studies (e.g. the closing of the mobil- 
ity gap [13]), the change of topology across this quantum 
phase transition has not been established theoretically, 
and the theoretical description for such a topological lo- 
calization transition is unavailable. Experimentally, how- 
ever, the disappearance of the topological FQHE with 
increasing disorder is common knowledge. 

In most studies of Anderson localization one increases 
the random disorder by fixing the strength of the impuri- 
ties (or the distribution of the strengths) and increasing 
the number of impurities. In this Letter, we adopt a 
different disorder model, which fixes the number of the 
impurities and increases the strength of each impurity 



potential. This construction provides a new path con- 
necting the gapped FQH state at weak disorder with the 
topologically trivial insulator at strong disorder. This 
enables us to investigate the phase transition(s) between 
these two topologically distinct insulators as well as to 
clearly distinguish the two phases. 

Although this new approach to disorder can be de- 
ployed in both lattice and 2D electron gas FQH systems, 
it is more natural to implement this in lattice systems for 
numerical studies. By increasing the potential strength 
of impurity sites in a lattice FQH model, this new ap- 
proach enables us to construct a complete phase diagram, 
which contains three different phases: a gapped topo- 
logical FQH state in the weak disorder limit, a trivial 
insulator phase at strong disorder, and an intermediate 
compressible gapless metallic liquid phase. To the best of 
our knowledge, such a quantum phase diagram has never 
been theoretically obtained in previous studies of FQHE. 
Even more remarkably, this construction allows us to ex- 
amine the nature of these phases and the phase transi- 
tions by measuring several different observables, includ- 
ing the energy spectrum, the entanglement spectrum, the 
Chern number, and the fidelity metric. Some of our find- 
ings (e.g., the lack of singularities at the topological trair- 
sition) offer important guiding principles for the general 
theoretical understanding of quantum topological phase 
transition. 

Model and phase diagram — The model we consider 
is a variation of the topological-fiatband model Q on a 
checkerboard lattice (Fig. [T]) with the Hamiltonian 



H = Hn 
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Here, the first three terms are the Hamiltonian of the 
flatband FQH model proposed in Ref. 0, 0] , and the last 
term describes the on-site impurity potential characteriz- 
ing the disorder strength, which is the tuning parameter 
for the quantum phase transitions. 

The first term (Hq) contains all the single-particle hop- 
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FIG. 1. (Color online) The checkerboard flatband model with 
charge impurities marked by the cones beneath the lattice. 

ping kinetic energy terms, 

Ho = tJ2 e'^'^ etc, + J2 t^A^J + t" E + h.c. 

Here, Ci [c]) is the fermion annihilation (creation) op- 
erator at site while and rep- 
resent the nearest-neighbor (NN), next-nearest-neighbor 
(NNN), and next-next- nearest- neighbor (NNNN) bonds, 
respectively. The U and V terms in Eq. ([1]) are the NN 
and NNN interactions with rii = c\ci being the par- 
ticle number operator. For the hopping strengths, we 
choose the same parameters as in Ref. [5[, i.e., t = 1, 
= ±7r/4, t' = 1/(2 + y/2) [f = -1/(2 + 72)) for the 
sohd (dashed) lines, and t" = l/(2 + 2-\/2). For a system 
with Nx X Ny unit cells (2 x iV^ x Ny sites), we choose 
NxNy = 3Ne, where Ne is the total number of fermions. 

The impurities describing the disorder are given by the 
last term in Eq. ([Ij. Here we select certain lattice sites 
where local potential traps are introduced to represent 
the impurities 

^impurity = ~ ^ Z^'*^'' (2) 

; 

where fii gives the strength of the potential trap at the 
impurity site / and the sum here is taken over the impu- 
rity sites only. We make a few non-essential assumptions 
to simplify the numerical study. We emphasize, however, 
that these assumptions do not affect our conclusions at 
all, and we have verified that relaxing them leads to the 
same qualitative results. We assume that the interaction 
between the impurity sites and the fermions is attractive 
{—fii < 0) and that all the impurity sites have the same 
fi. Due to this assumption (fj, being a constant), in order 
to achieve an insulating state with a unique ground state 
in the large /i limit, we need to require the number of im- 
purity sites Ni to coincide with the number of particles 
Ne- This condition Ni = is non-essential and is not 
needed if we relax the approximation of ji being a con- 
stant for all impurity sites. In the generic case where /i; 
varies from site to site, the large /i; limit is always a triv- 
ial insulating phase, regardless of the value of A^e, with 
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FIG. 2. (Color online) The phase diagram in the parameter 
space of impurity strength (^) and NNN interaction strength 
V for different system sizes {Ns = 2x3x4, 2x3x5 and 
2x3x6). Three phases are observed, the fractional-quantum- 
Hall phase (FQH), the compressible phase (metal), and the 
trivial insulator phase. The thick blue lines mark the averaged 
phase boundaries over different system sizes. Here we set 
f/ = 4 and iVe = Ns/&. 



all the particles pinned at the Ne sites with the largest 
/i;. To further simplify the numerics the impurity sites 
are chosen along one of the lines in our lattice as shown 
in Fig. m 

Using the exact diagonalization numerical method we 
establish the phase diagram of this model on finite lat- 
tices (Fig. [5]). The two extreme limits at ^ = and oo 
are easy to understand. For yU = 0, Eq. ([T]) recovers the 
interacting flatband model proposed in Ref. (without 
disorder). At 1/3 fiUing, this limit exhibits FQHE when 
suitable interactions are turned on 0, . For very large 
/i, the fermions are trapped at the impurity sites, result- 
ing in a topologically trivial insulator. For intermediate 
/i. a gapless compressible phase emerges as an interme- 
diate state, which leads to two phase transitions: (1) a 
topological transition between the topologically nontriv- 
ial FQH phase and the compressible phase, and (2) a 
typical metal-insulator transition between the compress- 
ible liquid and the trivial insulator. We note that for 
specific value of U and system sizes used in Fig. [21 while 
the transition from the FQH phase to the compressible 
metal can be driven either by increasing interaction {V) 
or disorder (/^), the transition to the trivial insulator can 
only happen for large disorder. 

Numerical measurements — The characteristics of 
these three phases can be seen from the energy spec- 
tra upon twisted boundary conditions, which are en- 
forced by requiring all the wave functions to satisfy 
T{Nj) I*) = e*"^ l^-), where j = a; or y and T{Nj) is 
the translation operator along the j direction. Here, we 
uniformly sample Ngx x Ngy points in the O^-Oy plane 
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FIG. 3. (Color online) Energy spectrum as a function of mag- 
netic flux (a) in the FQH phase with jj, — 0, (b) in the com- 
pressible phase with ^ = 2.5 and (c) in the trivial insulator 
phase with = 5. Here iVs = 2 x 3 x 6, iVe = 6, (7 = 4, 
and V = 1. And we choose Nex ~ Ney = 12. The blue 
dashed lines are the three lowest energies of momentum sec- 
tor Ky = ir/a with a being the lattice constant, and the black 
solid lines are the energies of the other momentum sectors. 



with 6,; 



2Tmgj /Ngj and ngj ~ 0,1, 



1. In 



Fig. |3] the energy spectra of the low energy states are 
plotted as a function of ngy x Ng^ + ng^ at different /i. 
At fj, — [Fig. mi^a)], a gapped state with three nearly- 
degenerate ground states is observed in agreement with 
Ref. , which is a characteristic feature of a topological 
FQH state. For fi = 2.5 [Fig.[3ljb)], there is no significant 
gap, indicating a gapless (compressible) phase. At large 
M [Fig- EIc)], we observe a gapped state with a unique 
ground state (i.e., a trivial non-topological insulator). 

To pin down the phase boundaries, we plot the excita- 
tion gaps averaged over twisted boundary conditions as a 
function of ^, Ai = Ei — Ei, where z = 1, 2, . . ., and Ei is 
the average energy of the state i. As shown in Fig. HJ^a), 
the FQH phase at small /i gives its way to the gapless 
phase at /i ~ 0.5. As /i is increased further, an excitation 
gap reopens at fj, ^ 3.6, where the compressible liquid 
turns into a trivial insulator. 

These conclusions are further supported by the entan- 
glement spectrum calculation, ^ = — 21n(pyi), which is 
obtained from the eigenvalues of the reduced density ma- 
trix PA = Trs/d with p being the density matrix of the 
ground state and Trs representing a partial trace over 
a part of the system. Here we split the system into two 
parts, A and B, in the particle space. Following the 
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FIG. 4. Numerical observables as a function of /i. (a) Aver- 
aged excitation gaps over twisted boundary conditions, (b) 
Particle entanglement spectrum. Here, we traced out 3 par- 
ticles, and the blue arrows mark the gaps in the entangle- 
ment spectrum. There are 75 states below the gap in the 
FQH state and this number becomes 10 in the trivial insula- 
tor phase, in good agreement with the counting of quasi-hole 
excitations, (c) Ground-state fidelity metric under periodic 
and twist boundary conditions, with 5^ = 0.001. In all these 
three panels, the red vertical dashed lines denote the phase 
boundaries and we set A^s = 2 x 3 x 5, A^e = 5, (7 = 4, and 
V = 1. 



method described in Ref. we find that at small /x, 
the entanglement spectrum contains a gap at ^ ~ 12, as 
shown in Fig. |4l(b), and the number of states below this 
gap agrees with the number of quasi-hole excitations in 
a FQH state. This gap disappears at fi 0.5 indicating 
a phase transition from a FQH state into a compressible 
state. At larger ^ (above ^ 3.6), a new gap emerges at 
^ ^ 7, and the number of states below this gap agrees 
with the number of excitations in a trivial insulator, sug- 
gesting the formation of a non-topological insulator. 

It is worthwhile to point out that due to the topological 
degeneracy, in order to get the correct entanglement spec- 
trum in the FQH phase, it is necessary to use the total 
density matrix pfqh = J2i which is summed 

over the three lowest states in the energy spectrum with 
i = 1, 2 and 3 Q. This is in direct contrast to the en- 
tanglement spectrum calculations in other phases where 
only the lowest state is used, p = |\I'i)(\I'i|. In or- 
der to present these two different calculations within the 
same plot, in Fig. SJb), we compute the entanglement 
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spectrum using an approximate density matrix, which is 
defined as 



1 ^ 



(3) 



Here Ei is the energy of the state i and the normahzation 
factor Af = X;Li exp[-/3(i;, - Ei)]. When 1/^ is much 
larger than the finite-size gap, but much smaller than the 
insulating gap, p recovers pfqh in the FQH phase and p 
in the trivial insulator phase within error bars. Any small 
difference between p and pfqh (or p) does not change the 
existence of the gap in the entanglement spectrum and 
the number of states below this gap. Therefore, p can be 
used to identify the nature of these phases for all values 
of p,. 

In addition, we also compute the Chern number in 
both the insulating phases, topological and trivial, us- 
ing flux insertion [18[. In the FQH phase, because the 
translational symmetry in our system is broken by impu- 
rities, the three nearly-degenerate ground states in our 
model can no longer be distinguished by different mo- 
menta. Therefore, it is not possible to observe the frac- 
tional Chern numbers C = 1/3 in this setup. Instead, the 
signature of such a FQH state is the total Chern num- 
ber being unity (similar to the case studied in Ref. jlTj). 
which is indeed what we observe here. In the trivial in- 
sulator phase, the Chern number is found to vanish as 
expected. This measurement directly probes the topo- 
logical structure of the two insulator phases (at small 
and large p), and demonstrates clearly the changing of 
the underlying topology as p increases. 

Fidelity metric and topological phase transition — Fi- 
nally, we present the ground-state fidelity metric g, which 
has been shown to be a sensitive indicator of quantum 
phase transitions in certain situations. The fi- 

delity metric g measures the change of the ground state 
wavefunction in response to a small change of the control 
parameter, defined as 
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2 l-F{p,Sp) 



[Spy 



(4) 



where F{p,Sp) = \{'i'{p)\'ii{p + Sp))\ is the overlap be- 
tween the two ground state wavefunctions at p and p+6p 
with 6p — > 0. The value of g remains small inside a 
quantum phase. However, at a quantum phase transition 
point, g diverges in the thermodynamic limit, reflecting 
the singularity associated with the quantum phase tran- 
sition. In a finite-size system, this divergence is usually 
regularized by the infrared cutoff resulting in a peak with 
finite height and width. The existence of such a peak is 
used as a signature of quantum phase transitions in nu- 
merical studies, and this method has been proven to be 
effective in detecting quantum phase transitions with the 
only known exception being the Kosterlitz-Thouless tran- 
sition in ID quantum systems, where the role of fidelity 



is still controversial [20|. In particular, it has been shown 
that this method can be used to detect topological transi- 
tions in the Kitaev model [2lj and in the integer quantum 
Hall effect [22j as well as in the topological transition be- 
tween a bosonic FQH phase and a charge-density-wave 
(CDW) state [H. 

As shown in Fig. IH^c), a peak at ^ ~ 3.6 marks the ex- 
pected phase transition between the compressible state 
and the trivial insulator phase. In addition, there are 
several smaller bumps at smaller values o{ p. In some 
delta peak may also appear (usually under high- 
symmetry boundary conditions with Ox and 0y being 
or tt). These bumps and delta peaks do not indicate any 
phase transitions. Instead, they come from the finite-size 
effect, and the existence as well as the heights and loca- 
tions of these bumps (delta peaks) are very sensitive to 
the boundary conditions {9x and 9y). In the FQH phase 
there are three degenerate ground states, whose ener- 
gies split due to finite-size effects. The values of these 
splittings (i.e., the finite-size gaps) are very sensitive to 
microscopic details such as boundary conditions and the 
value of p. For certain boundary conditions the finite- 
size gap between the lowest state and the first excited 
state may reach a minimum at some p, which will gener- 
ate a bump in g. In some special cases, these two states 
may cross at some p, and this level crossing will induce 



a delta peak in g as shown in Ref. [22|. Similar effects 
also arise in the compressible phase where the (finite-size) 
gap between the ground state and the first excited state 
is small. Ignoring these spurious bumps and delta peaks, 
we find no other feature or signature in g at the topolog- 
ical transition between the FQH state and the compress- 
ible liquid. The same is true for the excited states. We 
therefore conclude that our model reveals a new type of 
(topological) quantum phase transition which shows no 
singularity in q, in direct contrast to other topological 
transitions [H, [Hi, [ll] . 

Conclusion — We present a new approach to study the 
transition from a topological FQH state in weak disor- 
der to a topologically trivial insulator by increasing the 
strength of disorder. In our calculations, although the lo- 
cations of the impurity sites are not chosen randomly, the 
choice of their locations plays no qualitative role. There- 
fore, we expect that a topologically trivial random insu- 
lator obtained via Anderson localization is adiabatically 
connected to the topologically trivial insulator we find 
in our model at large p although our trivial insulating 
phase is probably more akin to a pinned quantum charge 
density wave phase. We emphasize that our approach 
can, in principle, be deployed in ordinary FQH systems 
in 2D electron systems. The only difference between the 
continuum case and the lattice system is that for the 
continuum model one needs to project all states to the 
lowest (or some other) Landau level to perform a numer- 
ical study, while in lattice models such a projection is 
unnecessary. In the presence of strong disorder and/or 
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strong interaction, such a projection is not a controlled 
approximation (for example, due to the strong Landau 
level coupling), which makes any numerical study in the 
continuum model less reliable. However, considering the 
strong similarities between the lattice model and the con- 
tinuum system, our qualitative phase diagram should in 
all likelihood apply to the continuum case also. 

This work is supported by DARPA-QuEST, Microsoft 
Q, and JQI-NSF-PFC. 



2, 389 (2011) 



[1] D. C. Tsui, H. L. Stormer, and A. C. Gossard, 

|Phys. Rev. Lettl 48, 1559 (1982) 
[2] H. L. Stormer, D. C. Tsui, and A. C. Gossard, 

|Rev. Mod. Piiysl n, S298 (1999) 
[3] C. Nayak, S. H. Simon, A. Stern, M. Freedman, and 

S. Das Sarma, |Rev. Mod. Phys?] 80, 1083 (2008) 

[4] E. Tang, J.-W. Mei, and X.-G. Wen, jPhys. Rev. Lett.| 

106, 236802 (2011) 
[5] K. Sun, Z. Gu, H. Katsura, and S. Das Sarma, 

|Phys. Rev. LetE] l06, 236803 (2011) 
[6] T. Neupert, L. Santos, C. Chamon, and C. Mudry, 

'Phys. Rev. Lett. 106, 236804 (2011) 

[7] D. Sheng, Z.-C. Gu, K. Sun, and L. Sheng. lNia~Comm.l 



[8] X.-L. Qi, .Phys. Rev. Lett.j 107, 126803 (2011) 
[9] N. Regnault and B. A. Bernevig, [Phys. Rev. X| l, 021014 
(2011) 

[10] S. A. Parameswaran, R. Roy, and S. L. Sondhi, 
arXiv:1106.4025 (unpublished) 

[11] G. Murthy and R. Shankar, arXiv: 1108.5501 (unpub- 
lished) 

[12] Y.-F. Wang, Z.-C . Gu, C.-D. Gong, and D. N. Sheng, 
|Phys. Rev. Lett] 107, 146803 (2011) 

[13] B. A. Bernevig and N. Regnault, arXiv:1110.4488 (un- 
published) 

[14] Y.-F. Wang, H. Yao, Z.-C. Gu, C.-D. Gong, and D. N. 

Sheng, arXiv:1110.4980 (unpublished) 
[15] Y.-L. Wu, B. A. Bernevig, and N. Regnault, 

arXiv:1111.1172 (un publishe d) 
[16] P. W. Anderson, Phy s. Rev.| 109, 1492 (1958) 
[17] D. N. Sheng, X. Wan, E. H. Rezayi, K. Yang, R. N. Bhatt, 

and F. D. M. Haldane, Phys. Rev. Lett. 90, 256802 

(2003) ^_ 

[18] Q. Niu, D. J. Thouless, and Y.-S. Wu, |Phys. Rev. B| 31, 

3372 (1985) 

[19] P. Zanardi and N. Paunkovic, ' Phys. Rev. ^ 74, 031123 
(2006) 

[20] S.-J. Gu, Int. J. Mod. Phys. B 24, 4371 (2010) 

[21] S. Yang, S .-J. Gu, C.-P. Sun, and H.-Q. Lin, 

|Phys. Rev. ^ 78, 012304 (2008) 
[22] C. N. Varne y, K. Sun, M. Rigol, and V. Galitski, 

|Phys. Rev. Bl 84, 241105 (2011) 



